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CHAPTER  I 
INTRODUCTION 

For  a  self-propelled  ship,  the  pressure  in  the  stern  region  is 
reduced  relative  to  that  for  a  towed  model  because  the  propeller  tends 
to  accelerate  the  water  flowing  into  it.  Due  to  this  action  of  the 
propeller,  the  forward  component  of  the  pressure  is  reduced,  increasing 
the  resistance  experienced  by  the  ship,  and  so  also  the  thrust  necessary 
to  propel  the  ship.  It  is  found  in  model  work  that  if  the  resistance  of 
a  hull  when  towed  at  the  speed  V  is  R,  the  thrust  necessary  to  propel 
the  ship  at  the  same  speed  V  is  greater  than  R  and  this  increase  in 
resistance  is  called  the  'augment  of  resistance'.  The  common  practice 
is  to  look  upon  this  increase  in  resistance  as  a  deduction  from  the 
thrust  provided  by  the  propeller,  so  that,  although  the  propeller  pro¬ 
vides  a  thrust  of  T, tons,  only  R  tons  are  available  to  overcome  the 
resistance.  This  loss  of  thrust  (T-R) ,  when  expressed  as  a  fraction 
of  the  thrust  T  is  called  "the  thrust  deduction  fraction,  t"  where 

t  =  (T-R) /T  =  1  -  R/T. 

In  a  recent  investigation  by  Huang  et  al  [1},  one  of  the  conclusions 
was  that  "Potential-Flow  Propeller-Hull  interaction  methods  predict 
thrust  deduction  and  propeller-induced  pressure  distribution  very  well." 
The  potential- flow  methods  used  in  [1]  are  based  on  the  well-known 
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and  widely  used  computer  program  of  Hess  and  Smith  [2] .  It  was  found 
that  the  pressure  drag,  for  irrotational  flow  without  propeller,  deviates 
from  the  theoretical  value  of  zero  by  an  amount  of  the  same  order  of 
magnitude  as  the  viscous  pressure  drag,  Ryp.  This  drag  error  is  a 
measure  of  the  accuracy  of  that  computer  program.  Nevertheless,  the  H.-S. 
method  was  also  used  to  compute  the  pressure  distribution  and  pressure 
drag  with  propeller,  and  the  difference  between  the  so-computed  pressure 
drags  was  assumed  to  give  the  thrust  deduction.  Also,  it  is  known  that 
the  thrust  deduction  is  of  the  same  order  of  magnitude  as  the  viscous 
pressure  drag,  R^p,  and  hence  this  method  of  computing  thrust  deduction 
could  lead  to  large  percentage  errors.  The  validation  of  theoretical 
procedures  of  [1]  requires  an  irrotational-flow  calculation  method  more 
accurate  by  at  least  an  order  of  magnitude  than  that  of  H.-S.  to  resolve 
the  small  difference  between  propeller  thrust  and  viscous  pressure  resis¬ 
tance. 

In  order  to  determine  the  thrust  deduction  when  the  boundary  layer 
and  near  wake  are  taken  into  account,  it  is  necessary  to  calculate  the 
flow  characteristics  under  these  conditions,  with  and  without  a  flow- 
producing  mechanism  simulating  the  presence  of  a  propeller.  For  accom¬ 
plishing  this,  it  was  proposed  to  make  use  of  a  computer  program  for  cal¬ 
culating  the  thick  boundary  layer  by  a  differential  method.  Such  a  program, 
developed  by  Lee  and  Patel  [3]  for  axisymmetric  flow  past  a  body  of  revo¬ 
lution,  was  available,  but  some  modifications  were  required  for  the  present 
purpose.  Regarding  the  irrotational  flow  calculations,  L.  Landweber  and 
D.H.  Choi  have  developed  a  program  which  calculates  the  velocity  distri¬ 
bution  on  the  surface  of  the  body  to  the  required  accuracy  for  drag 
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calculations.  Thus  the  means  for  calculating  thrust  deduction  for 
a  body  of  revolution  by  irrotational-flow  methods  and  from  the  boundary- 
layer  equations  were  at  hand.  Hence  this  investigation  to  confirm  the 
validity  of  the  irrotational-flow  procedure  was  undertaken. 

It  is  proposed,  first,  to  solve  the  thick  boundary- layer  equations 
for  uniform  flow  past  the  body  by  a  differential  method  [3]  to  find  the 
drag  of  the  body.  Next  the  thick  boundary- layer  equations  will  be 
solved  to  find  the  drag  with  a  mechanism  representing  the  effect  of  the 
propeller.  The  outer  irrotational-flow  equations  will  be  rederived  in 
Chapter  III  with  an  additional  singularity  representing  the  effect  of 
the  propeller. 
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CHAPTER  II 

SELECTION  OF  FORM  AND  SINGULARITY 


The  thick-boundary -layer  equations,  when  used  to  compute  the  boun¬ 
dary-layer  characteristics  for  the  modified  spheroid  [3],  broke  down 
at  a  point  where  there  is  discontinuity  in  the  surface  curvature.  In 
order  to  resolve  this  problem,  a  more  detailed  analysis  of  the  Reynolds 
equations  is  required.  Hence  a  streamlined  body  of  revolution  with  con¬ 
tinuous  slope  and  curvature  was  selected.  The  shape  adopted  is  the  EPH 
form,  so-called  because  of  its  elliptical  nose,  parabolic  after  body,  and 
hyperbolic  tail.  This  was  proposed  by  L.  Landweber,  and  adopted  in  about 
1945,  as  a  cavitation  resistant,  low-drag  strut  for  the  U.S.  Navy.  The 
equation  of  the  form,  in  cylindrical  coordinates  (x,r),  is 


with  a/b  =  5.5  and  L/2b  =  6.3055.  Here  a  is  the  length  of  the  ellipti¬ 
cal  nose,  b  is  half  the  maximum  diameter  of  the  body  and  L  is  its  over¬ 
all  length.  See  fig.  1.1. 


1.  Selection  of  Mechanism  to 
Represent  a  Propeller 
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It  is  common  practice  to  use  a  sink  disc  behind  the  tail  in  order 
to  represent  the  propeller.  For  the  present  purpose,  however,  a  sink 
was  chosen  in  order  to  simplify  calculations,  since  our  purpose  is  to 
investigate  a  principle,  rather  than  to  perform  a  design  calculation. 

The  position  and  strength  of  the  sink  can  be  selected  to  match  the 
thrust  deduction  of  the  propeller. 

The  position  of  the  sink  on  the  axis  behind  the  tail  is  chosen  arbi¬ 
trarily  at  half  the  maximum  beam  from  the  tail.  Having  fixed  the  posi¬ 
tion  of  the  sink,  it  now  remains  to  determine  the  strength  of  the  si'.ik 
that  gives  the  required  pressure  drag  in  potential  flow,  i.e.,  the  cal¬ 
culated  thrust  deduction.  For  this,  the  drag  coefficient  Cp  is  calcula¬ 
ted  by  integration  of  the  pressure  on  the  body  surface  and  an  expres¬ 
sion  for  CQ  in  terms  of  the  strength  M  of  the  sink  is  obtained. 

The  pressure  drag  D  in  irrotational  flow,  in  which  there  is  a  uni¬ 
form  stream  of  unit  velocity  past  the  body  and  a  sink  of  strength  M  on 
the  axis  behind  the  tail,  is  given  by 


D  =  /  p  2irr 
0 


dr 

dx 


dx 


(2.1) 


where  p  is  the  pressure  on  the  surface  of  the  body,  and  r  is  its  ordi¬ 
nate.  This  may  also  be  written  as 


D 


0 


(P-P«l 


d(r2) 

dx 


dx 


(2.2) 


j  ar/iij  .zit'hf*. 


idierep^  denotes  the  pressure  at  infinity.  If  we  define 
2  2 

C  =  20/CpU^  A)  and  C  =  2(p-peo)/pUoo  where  is  the  free-stream  velo- 
1/  p 

2 

city,  A  is  the  maximum  cross-sectional  area  ird  /4,  and  d  is  the  maximum 
diameter,  then  equation  (2.2)  can  be  written  as 


d(r2) 


c  =  f  C  5LL1 
CD  d2^0Cp  dx 


(2.3) 


We  also  have,  from  the  Bernoulli  equation, 


C  *  1  -  —r 
P  n2 


(2.4) 


or,  expressing  U  in  the  form 


U  =  T(x)  cos  a(x) 


(2.5) 


where  U  is  the  tangential  velocity  and  a  the  angle  of  the  tangent  with 

2 

the  x-axis  at  a  point  on  the  body,  and  putting  r  =  f(x),  we  obtain 


U  =  T(x)  [4f/(4f+f »2) ] 1/2 


(2.6) 


Using  (2.6)  and  (2.4),  (2.3)  may  be  written  as 


cn  =  ‘  ~TJ  A2W  [ff 7(4f  +  f’2)]dx 

uV  0 

00 


(2.7) 


The  above  equation, which  is  valid  for  axisymmetric,  irrotational  flow 

about  a  body  of  revolution  immersed  in  a  uniform  stream,  should  give  the 
value  zero. 
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We  shall  now  derive  a  form  of  CQ  when  an  axial  downstream  sink  of 
strength  M  is  present.  Since  the  boundary  condition  and  the  Laplace  equa¬ 
tion  are  homogeneous  and  linear  superposition  of  solutions  is  valid, 
r (x)  may  be  assumed  to  be  of  the  form 

T(x)  =  UJ1  +  Mr2  (2.8) 

_2 

^ere  F^  is  nondimensional ,  but  has  the  dimensions  L  .  Substituting 
for  r(x)  in  (2.7)  and  noting  that,  with  the  uniform  stream  alone,  the  drag 
is  zero,  we  now  obtain 

Cn  =  T2/  (2MU/  r  +  M2rb[ff'/(4f+f'2)]dX  (2.9) 

D  ua  0  1  2 

00 

Hence  the  drag  coefficient  Cp  can  be  written  as 

CD  =  AM’  +  BM'2  (2.10) 

M 

where  M*  =  - and  A  and  B  are  the  nondimensional  coefficients 

U  d^ 

00 

L  ? 

A  =  -  32  /  r  r„[ff'/(4f+f'  )]dx 

0  1 

(2.11) 

B  =  -  16  d2  /  r2  [ff ’/ (4f+f ,2)]dx 
0 

The  coefficients  A  and  B  of  (2.9)  were  calculated  from  (2.11)  by  deter¬ 
mining  r2  and  r2  from  the  aforementioned  program  of  Landweber  and  Choi. 

The  results  obtained  are 
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A  =  -  0.98195474 
B  =  0.7907678S 

Thus  CQ  is  given  by  a  quadratic  equation  in  M' ,  (2.10). 

The  approximate  drag  coefficient  of  a  low  drag  body  at  a  Reynolds 
number  of  the  order  of  106,  is  about  0.04.  We  also  know  that  the  viscous 
pressure  drag  is  approximately  0.1  times  that  of  the  total  drag  and  is 
also  approximately  equal  to  the  thrust  deduction.  Hence  a  drag  coeffi¬ 
cient  CD  of  0.004  was  used  in  the  present  calculations.  Solving  the 
quadratic  equation  for  Ml  we  obtain  the  value  of  M*  =  -0.004  for  the 
sink;  the  other  value  for  M'is  neglected  since  it  is  positive. 
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CHAPTER  III 
IRROTATIONAL  FLOW 


The  irrotational  flow  outside  the  boundary  layer  and  wake  can  be 
derived  as  a  solution  of  an  exterior  Neumann  problem  of  potential  theory. 
As  is  well  known,  this  type  of  problem  can  be  solved  by  setting  up  an 
integral  equation  of  the  second  kind  in  a  surface  distribution  of  singu¬ 
larities.  In  the  present  problem,  an  integral  equation  of  the  first  kind 
is  considered,  for  three  reasons.  Firstly,  the  calculation  is  relatively 
simple  and  less  time-consuming  in  comparison  with  methods  based  on  an 
integral  equation  of  the  second  kind.  Secondly,  the  boundary  conditions 
can  be  prescribed  in  the  form  of  the  stream  function  on  the  'fixed 
matching  boundary'  and  the  solution  of  the  integral  equation  yields  dir¬ 
ectly  the  velocity  distribution  on  the  boundary  which  is  required  as  an 
input  for  the  boundary- layer  and  wake  calculation.  This  avoids  addi¬ 
tional  operations,  such  as  numerical  integration  or  differentiation,  on 
the  solution  of  the  integral  equation  for  the  subsequent  boundary- layer 
calculation  to  obtain  the  necessary  boundary  conditions  for  the  poten¬ 
tial-flow  calculations.  Thirdly,  the  high  accuracy  given  by  the  inte¬ 
gral  equation  of  the  second  kind  is  not  required  for  boundary- layer 
calculations. 

For  the  present  case,  with  a  sink  on  the  axis  behind  the  tail,  the 
integral  equation  has  to  be  rederived.  For  this  purpose  let  us  consider 
the  general  case  of  a  body  introduced  into  a  flow  field  of  undisturbed 
velocity  potential  4>0* 
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1.  Irrotational  Flow 
Outside  the  Closed  Body 

See  (Fig.  3.1)  The  total  velocity  potential  $  can  be  written  as 


+  4 


(3.1) 


where,  <J>q  denotes  the  velocity  potential  of  the  ambient  flow  field  and 
<t>  the  disturbance  potential  due  to  the  introduction  of  the  body.  The 
boundary  condition  on  the  surface  S  of  the  body  is 


3$  _  9*o  ,  34>  _  0 
3n  3n  3n 


(3.2) 


where  n  denotes  the  distance  normal  to  S  from  a  point  on  S.  Let 
3  2  3 

4*^  *  -(x-t)/R,  =  r  /R  denote  the  potential  and  stream  functions  of 

a  doublet  of  unit  strength  oriented  in  the  positive  x-direction  at  a 
point  (t,0)  on  the  axis,  and 

R  =  [(x-t)2  +  r2(x)]1/2  (3,3) 

2 

The  equation  of  the  body  is  r  =  f(x),  0  x  <_  L  in  cylindrical  coordi¬ 
nates  (x,r).  Let  U  (x)  be  the  tangential  velocity  at  a  point  on  the 
body  surface,  s  the  arc  length  measured  from  the  nose  N  of  the  body, 
and  l  the  length  of  arc  from  the  nose  to  the  tail  along  a  meridian.  By 
Green’s  reciprocal  theorem,  we  have 


(3.4) 
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where  dS  is  the  surface  area  element.  By  (3.2)  and  (3.4)  we  have 


l  3<j>  , 
/  r*  d 
0 


3n 


ds 


i  3(f) 

0 
o 


"d  3n 


ds 


(3.5) 


where  ds  is  the  element  of  arc  length.  Since 


3n 


^d 

3s  ’ 


d<P  d\p 


3n 


3s 


(3.6) 


where  \pQ  denotes  the  stream  function  corresponding  to  $o,  (3.5)  then 
becomes 


Z  dtp 


I 


ds 


ds  = 


Z 

-/ 


dtp 

»d  te¬ 


ds 


By  integration  by  parts,  we  obtain 


M- 


l  l  4 

'  I  *d  M  ds  *  '  /  ♦ 


3*, 


o  o 


d  3s 


ds 


Since  r  =  0  at  s  =  0  and  s  5  H,  s  0  at  both  ends,  and,  by  substitu- 

3<j>  3*  d 

ting  a? =  I!  -  jr  * u  -  t r-  "e  oUain 


l 

/U  3 

o  R 


2.  .  I  3<p  dtp 

^  da  =  /t*d^*  *aIT)da 


(3.7) 


Equation  (3.7)  is  in  the  form  of  an  integral  equation  of  the  first  kind 
for  Ug  since  the  right  member  is  a  known  function. 
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It  can,  however,  be  shown  that  [tj^  d<j>Q  +  <fyo]  is  an  exact  dif¬ 
ferential  (see  Appendix) .  Consequently  the  integral  of  the  right  hand 
side  of  equation  (3.7)  is  independent  of  the  path  of  integration  between 
nose  and  tail  points  of  the  body.  Hence  the  integral  may  be  evaluated 
along  the  indicated  path  in  Fig.  3.2,  along  the  axis  from  nose  N  to  the 
tail  T,  with  the  doublet  at  (t,0)  avoided  by  means  of  a  semi-circle  of 
very  small  radius  R. 

Along  the  axis,  we  have  =  0,  and  since  \pQ  is  an  axisymmetric 
stream  function,  tpQ(x,0)  is  a  constant  along  the  axis.  Hence,  we  need  to 
evaluate  the  right  member  of  (3.7)  only  along  the  small  semicircle.  Then 
the  integral  becomes 


J  [* 


3<t 


*0, 


+  A  —Tl  (JS 

d  3s  fd  3s  J 


ir  .  2  3<J> 

t  .sin  6  0 

j  1  R  36 


cos6 
r.2  30 


]  de 


where  use  has  been  made  of  the  relations  r  =  R  sin  0;  x  -  t  =  R  cos  0 
and  R  is  very  small.  We  have  for  and  -rr~, 

00  do 

3^0  3*o  n  3*o 

W  =  -  R  sine  35T  +  R  cose  TT 

(3.8) 


dip  dip  dip 

W’  -  R  sin9H£*  R  cose  3T 


3*  S*o 

=  R  sine  r  - —  +  R  cos0  r  - 

3r  3x 


5V  .  ,8V 


We  now  expand  and 


in  Taylor  series  about  the 


point  (t,o) : 
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3  4> 

'■Sx  '  (x,r) 


34>0 

f — — ) 

'•ar  J(x,r) 


2  2 

3<j>  3  4>  3  <p 

^lx~\,0  +  Cx-1:)  (  2)  t>0  +  r  C3x3r)tj0  +  - 

oX 

3$  3  4>  3^(j> 

..  O.  .  w  O.  ,  CK 

C3r  ?t,0  +  (x_t)  (3x3r5t,0  +  r  (  2  ^t,0  + 

or 


(3.9) 


Using  the  fact  that  (^)  =  o  by  symmetry  in  (3.9)  and  neglec¬ 
ting  higher  order  terms,  we  obtain  from  (3.8) 


3<f>  3<j> 

jr  ”  - R  si"9  fir’t.o 


3»|>  2  3<|> 

aT  ’  R  sln9  cose  ‘ir5  t.o 


(3.10) 


Hence  the  integral  becomes 


X.  L  cJffi 

/  U  ~  ds  =  2  (-t-^1  , 
i  s  d3  l3x  Jt ,( 


(3.11) 


If  4>o  =  U^x  where  U^  =  1,  equation  (3.11)  reduces  to  the  well-known  inte¬ 
gral  equation  for  Ug  on  the  surface  of  an  axi symmetric  body  placed  in  a 
uniform  stream  of  unit  velocity. 


a  2 

f  U  ~  ds  =  2 
o  s  R 


2.  Irrotational  Flow  Past  the  Semi-Infinite 
Body  Enclosing  the  Boundary  Layer  and  Wake 

Spis  a  fixed  matching  boundary,  (FMB)  enclosing  the  boundary  layer  and 

wake,  s  is  the  arc  length  along  a  meridian  measured  from  the  nose  N,  n  is  the 

distance  normal  to  Spfrom  a  point  on  S^and  V  the  value  of  the  stream 


14 


function  corresponding  to  the  total  velocity  potential  4,  is  pre¬ 
scribed  along  s.  See  fig,  [3,3  ].  The  boundary  condition  on  the  sur¬ 
face  is  given  by 


3$  _  8*o  t  3<j> 
3n  3n  3n 


or 


(3.12) 


3 V  _  *o  +  3 £ 
3s  3s  3s 


Applying  Green's  reciprocal  theorem  to  potentials  and  (j>  gives 


3<j> 


/  *  aif*  ds  =  /  ds 


d  3n 


(3.13) 


where  dS  is  the  surface  area  element.  From  (3.12)  and  (3.13)  we  obtain 


3<j> 


co  8(b 

3$  ^ 


/  r  ♦  ds  =  /  r  ^  ~  a^2-]  ds 


dl3n  3n 


(3.14) 


3$ 


3y 


Using  (3.6)  and  the  relation  r  — =-■» —  in  (3.14),  we  find 

oil  oS 


00  ”  aw  00  3,L 

/  ♦  zT  ds  =  /  ^  h  ds  -  /  —  ds 


fd  3s 


Kd  3s 


or,  after  integration  by  parts,  the  above  equation  may  be  written  as 


i  H*dds  = 


00  oo  00 


'd  3s 


o 


15 

3<j>  o<J) 

Since  =  - -r— -  =  U  -  t — »  r  =  0  at  s  =  0,  and  # ,  as  R-*»  is  zero 

oS  oS  oS  oS  d 

this  becomes 

«  2  “  .  00 

/  »  ^ ds  =  -  H  i ds  *  /  [*h  sr  *  ^  jri ds  <3-ls> 


R 


The  integral  /  [<j>^  ^  -^— ]  ds  is  independent  of  the  path  of  inte- 

3<f> 

gration  and,  as  was  shorn  in  the  previous  section,  is  equal  to  2(- — )  . 

uX  t  j  u 

Hence,  we  obtain 


r  „  r2  .  r  A  3» ,  ,/v 

I  u  I?  ds  m  /  *d  sTds  *  2(to_) t,o 

o  R  o 


or 


(3.16) 


00  2  “  3(j)  ,  3<J> 

I  U  3ds  =  / 

o  R  o 


Since  'i'(O)  =0  and  -*  0  as  R  -**>. 

3.  Velocity  at  Matching  Boundary  in  Presence  of 
an  External  Sink 

The  interaction  of  the  boundary  layer  and  wake  with  irrotational 
flow  is  solved  by  an  iteration  procedure  [3] .  It  is  now  necessary  to  take 
into  account  the  effect  of  an  external  sink  on  these  calculations,  since 
a  sink  is  chosen  for  the  present  "thrust  deduction"  calculation  (Ch.  II). 
The  normal  velocity  to  the  FMB  and  the  value  of  the  stream  function  on 
the  FMB  due  to  an  external  sink  and  its  image  system  within  the  body 
are  required.  The  normal  velocity  and  stream  functions  at  the  FMB, 
calculated  by  the  interaction  procedure  for  the  flow  without  the  ex¬ 
ternal  sink,  will  now  be  modified  by  the  contributions  of  the  sink 
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and  its  associated  disturbance  potential  in  the  potential  flow.  These 
updated  values  of  normal  velocity  and  stream  function  on  the  FMB  are 
used  to  calculate  the  external  potential  flow  (outside  BLW)  once  more, 
and  the  resulting  pressure  distribution  is  used  to  calculate  the  BLW. 


a)  Method  1:  Using  Green's  third  formula.  The  tangential  velocity 
U  on  the  surface  of  the  body  is  obtained  by  solving  a  Fredholm  integral 
equation  of  first  kind,  in  the  presence  of  an  external  axial  sink.  The 
disturbance  velocity  potential  on  the  surface  of  the  body  can  then  be 
written  as 


(x,r  )  =  J  U  ds 


(3.17) 


M 

*'0v'v,ibJ  ”  J  u  as  -  pp 
^  O  S 

where  here,  and  in  the  following,  Q(S,n>C)  denotes  a  point  on  the  body 

2  2  21/2 

surface,  M  the  sink  strength,  Rg  =  [(5-xs)  +  n  +  C  ]  is  the  distance 

from  the  position  S(xs,0)  of  the  sink  to  a  point  Q  on  the  body. 

With  this  value  of  the  disturbance  potential  on  the  surface.  Green's 
third  formula  can  be  used  to  calculate  the  disturbance  potential  <|>  at 
any  point  P(x,y,0)  in  the  field.  We  have,  by  Green's  third  formula. 


= 


/  [*r 


i 


p  4n  '  Q  anQ  Rpq  RpQ  3nQJ  Q 
where  RpQ  =  [(x-C)2  +  (y-n)2  +  ;2]1/2  4 


]  dSQ  (3.18) 

is  the  distance  between  the  points 


P  and  Q.  [See  fig  (3.4)]. 

The  boundary  condition  on  the  body  surface  is  that  the  normal  velo¬ 
city  is  zero. 


3<|>  3<j> 

*  +  -z —  =  0  on  S,  <ji 

3nQ  3nQ  ‘‘s 
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or 


9<J> 


3n, 


Q  _ 


=  -  M 


3  1_ 

3„  K 

nQ  s 


(3.19) 


We  have 


and 


4  k  ~~  if1  f'(v6HQ  * 11  "q  * 

^  s 


anQ  Rpp  =  -3  [(X"5HQ  +  (y"n>Q  "  ?nQ3 

where  £p,  np  and  np  are  the  direction  cosines  of  the  normal  to  the  body 
at  the  point  Q  and  are  given  by 


=  -  =  •—  and  D, 


D, 


Q  V  Q  °Q 


=  A'2  + 


4f 


With  the  help  of  relations  for  i,  m,  n  and  noting  that  n  =  rp  cos  0p 
r;  =  Tq  sin  0p  and  y  =  rR  when  0p  =  0  in  cylindrical  coordinates,  we 
can  rewrite  the  above  equations  as, 


3  1 


3nQ  Rs  DnR 

x  Q  s 


H  ^xs  "  ^£'  +  2tq} 


(3.20) 


and 


-  -  77T  tc*  -  o*’  -  2Vq  cos  9q  *  2rQ2i 


Q  ^  Vpq 


(3.21) 
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Consider  a  source  situated  at  the  point  P.  Then  by  Gauss  flux 
theorem,  we  have 


dS„  =  0 


PQ 


(3.22) 


Equation  (3.22)  is  used  to  remove  the  peak  in  the  integral  of 
equation  (3.18)  by  writing  it  in  the  form 


1_ 
4  IT 


^  [(vv5  3*V  R 


3.  1 

Q  “PQ 


1  , 
RPQ8nQ  Q 


(3.23) 


where  <J>pi  is  the  potential  at  the  point  P',  the  point  on  the  body  such 
that  PP'  is  normal  to  the  body.  When  Rp^  is  small,  (<f>  q'^pf)  is 
also  small. 

Equation  (3.23)  may  be  written  as 


^P  "  4rr  F  {(4>Q  '  V5  L  sL  R, 


d9. 


0  °  0  PQ 


3*0  ^  1 

-a r  /r-d0o}ds  (3-24) 

>  0  RPQ  Q 


where  P  is  the  distance  from  nose  to  tail  measured  along  a  meridian 
and  ds  is  an  element  of  arc  length.  We  have,  from  equation  (3.21), 


2ir 


3  1 


S0  *"q  S  de  ” 


i  ,  2u  ,  2rrtrn  2 it 

-  jlp-  [(x-C)f'  +  2r  2]  j  ~  d9  ♦  /  S2|®  d 9  (3.25) 

2  Q  0RPQ  **  0RPQ 


We  now  need  to  evaluate  the  integrals 


I,  = 


2rr  , 

/  -S-de,  I. 

0RPQ 


2n 


0  R, 


i. 


-  I 


2tt 


d  9 


PO 


0  RPQ 


We  can  write 


%  ■  <X<  '2  *  rp  '  *  Tq  -  2rprQ  C0S 


where 


=  (x-g)2  +  (rp+rQ)2  -  4rPrQ  cos2  6/2 


R'2  [1-k2  cos2  -£] 


R'2  =  (x-g) 2  +  (rp+rQ)2,  k2  =  4rprQ/R'2 


Substituting  3=  —  (tt-0)  in  the  above  integrals,  we  obtain 

I  =  f_L.  d6  =  _i_  Y2 _ _  =  _4_  EOO 

1  0  rJ  R'3  0  fl-k2sin2Sl  3/2  R'3  1-k2 


2ir 

t  -  r  cosq 

a2  ■*  3 

0  rpq 


.  tt/2  ,  .  .2 

de  =  -  Ar  J 


“T  ~  "I— -P-l/2-  dS 

R,-i  0  [1-k  sing] 57 ^ 


=  1t^-*m4  =  -Vtkw 

R'-5  K  l-k^  R'V 

2tt  ,  .  tt/2  _ 1 _ 

b  ■  l  [l-*2sinV/2 


where  K(k)  and  E(k)  are  the  elliptic  integrals  of  first  and  second  kind 
respectively,  he  have 


?  4rprn 

1-k  =  1 - ^ 

(x-g)‘ 


2 

=  !*l 

R-2 


(3.26) 


where 


rpq  =  (x_5)2  +  (VrQ)2  and 


2  2  "> 
2-k  =  _1 _  2_  R’  R'" 

k2(l-k2)  1-k2  +  k2  Rp^  2rprQ 


(3.27) 


Substituting  (3.26)  and  (3.27)  in  I  ,1  and  I  and  rewriting  (3.24) 
with  the  help  of  and  (3.21),  we  obtain,  after  collecting 

terms, 

[2rQ(rp-rQ)-(X-Of'(0] 
ro 

~  u:  (3.28) 

Q 

We  have  ds  =  dxsec«  =  ^—7  where  a  is  the  angle  made  by  the 

Q 

tangent  with  the  x-axis.  4l  can  now  be  written  as 


*p  =  Pi iP  *(E-K)  +  irL*3l2rQ(VroHx^f[?)] 

1  rpq 


^fr[(V5)P*  2rQ2)jdx 


(3.29) 


where  the  body  extends  between  x  =  -1  and  1. 

Equation  (3.29)  is  used  to  calculate  the  disturbance  potential  at 
five  closely-spaced  points,  two  on  either  side  of  the  fixed  point 
where  the  normal  velocity  to  FMB  is  desired,  and  one  coinciding  with 
the  fixed  point.  The  five  points  are  taken  along  the  normal  to  FMB 
at  this  point.  The  normal  derivative  of  the  potential  is  obtained  by 
fitting  a  Lagrange  interpolation  polynomial  through  these  five  points. 
The  normal  derivative  of  the  disturbance  potential  yields  the  normal 
disturbance  velocity  V  *n.  To  obtain  the  total  normal  velocity,  that 
due  to  the  sink  is  to  be  added  to  the  above  value  of  V  *n. 

The  velocity  at  a  point  P  due  to  a  sink  of  strength  M  placed 
at  a  point  (xs>0,0)  is  given  by 
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Having  calculated  the  normal  velocity  V  the  stream  function 
iJj  at  any  point  can  be  calculated  by  integrating  V  along  a  meridian  from 
the  nose 

s 

i|»  =  -  /  r  V  ds  (3.32) 

0 

The  tangential  velocity  U  on  the  body  surface,  used  as  an  input  to 
Green's  third  formula  in  calculating  the  velocity  potentials  and  hence 
the  normal  velocities  at  the  FMB,  is  obtained  by  solving  an  integral 
equation  of  first  kind  in  U.  This  velocity  is  accurate  upto  four  signi¬ 
ficant  digits  only.  Due  to  the  nature  of  integral  equations  of  the  first 
kind,  greater  accuracy  cannot  be  expected  from  this  solution.  When  this 
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solution  was  used  in  Green's  third  formula  to  calculate  velocity  potential, 
it  failed  to  give  a  good  estimate  of  the  normal  velocity  at  FMB.  This  is 
so  because  the  calculation  of  normal  velocity  involves  diffe  ^nces  in  velo¬ 
city  potentials  of  nearly  equal  magnitudes.  This  became  apparent  when  the 
computer  program  was  checked  with  the  exact  solution  for  a  sphere.  The 
use  of  Green's  third  formula  demands  inputs  that  are  more  accurate  than 
that  given  by  the  solution  of  an  integral  equation  of  the  first  kind.  As 
a  result,  it  was  decided  to  try  a  different  approach  to  obtain  both  normal 
and  tangential  velocities  due  to  the  sink  and  its  image  systems  at  FMB,  as 
presented  in  the  next  section. 

b)  Method  2:  Using  Taylor  series.  By  this  method,  we  can  not  only 
estimate  the  normal  velocity  at  FMB  due  to  the  sink  and  its  image  system, 
but  we  can  also  estimate  the  tengential  velocity  at  FMB  due  to  the  sink 
and  its  image  system.  This  eliminates  the  need  to  solve  an  integral  equa¬ 
tion  for  flow  past  the  semi-infinite  body  including  the  effect  of  the  sink. 
Instead,  the  calculated  normal  and  tangential  velocities  due  to  the  sink 
and  its  image  system  within  the  body  are  added  vectorial ly  to  those  ob¬ 
tained  by  solving  an  integral  equation  (3.16)  for  uniform  flow  past  the 
semi-infinite  body. 

The  resulting  pressure  distribution  is  used  to  calculate  the  turbu¬ 
lent  boundary- layer  once  more,  thus  incorporating  the  effect  of  the  sink 
on  the  calculation  procedure. 

Equation  (3.11)  was  solved  to  obtain  the  surface  distribution  of 
tangential  velocity  u(s,0)  on  the  closed  body.  Since  the  FMB  is  not  very 
far  from  the  body,  it  is  assumed  that  the  normal  to  the  body  is  also 
normal  to  the  FMB.  Now  a  procedure  for  estimating  normal  and  tangential 
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velocities  at  FMB  due  to  the  sink  and  its  image  system  will  be  presented. 


By  Taylor  series,  we  can  write  for  U  and  V 

2 

U(s,n)  =  U(s,0)  +  n  Ug  (s,0)  +  ^  USSCS»°)  +  .  (3.33) 

V(s,n)  =  V(s ,0)  +  n  Vs  (s,0)  +  ^  vss(s»°)  +  •••  (3.34) 


By  noting  that  the  normal  velocity  on  the  body  surface  is  zero  and 
neglecting  higher-order  terms,  we  can  write 


2 

U(s,n)  =  U(s,0)  +  n  Us(s,0)  +  ^  Uss(s,0) 
V(s,n)  =  n  Vs(s,0)  +  — ■  Vss(s,0) 

We  have  the  equation  of  continuity 


(3.35) 

(3.36) 


a?  ™  *  h  (hirV)  - 0  (3-37> 

Since  the  flow  is  irrotational ,  we  have 


§  -  h  (hiu)  -  0  ‘3-38> 

where  r  =  rfe(s)  +  n  co£  r(s) ,  hj  =  1  +  nK(s)  is  the  longitudinal  metric 
coefficient,  K(s) ,  a(s)  and  r^(s)  are  the  longitudinal  curvature  of  the 
body,  the  angle  of  the  tangent  to  the  body  with  the  x-axis  and  the  body 
ordinate  at  the  point  (s,0)  on  the  body,  respectively.  In  the  following 
equations,  the  notation  U(s,n),  V(s,n)  etc.,  has  been  omitted  and  U,V, 
etc.,  will  be  used  for  simplicity.  Eqns.  (3.37)  and  (3.38)  ®ay  be 


written  as 
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W?*1!*' #  'Ji* *■  >S^»  'ft**-***-'  -s-rr-^j^ .  __ 
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b  Us0  *  Vina  *  rbVno  =  0 

(3.45) 

'nO  *  V  "  0 

(.'.to) 

*  U  A  +  r.V  _  +  2U  n  sina  =  0 
b  ssO  b  snO  sO 

:.  U  A  +  U  A  cosa+  U  A  sina+  KU  sina  +  r,V  _  + 
b  snO  sO  nO  o  b  nnO 

(3.47) 

2V  A  (Krn  +  cosa)  =  0 
nO  0 

(3.48) 

V  .  +  KU  n  +  K'U  =  0 
snO  sO  0 

(3.49) 

V  n-U  n  -  KU  =  0 
snO  nnO  nO 

(3.50) 

By  (3.45),  we  have 

V  (s,0)  =  vrn  =  -  UA  ^2.  -  u 
n  n0  0  r,  sO 

b 

(3.51) 

From  (3.46) , 

Un(s,0)  =  Un()  =  -  UoK 

(3.52) 

From  equations  (3.48),  (3.45),  (3.46)  and  (3.4g),  we  obtain 


V(S,0)  “  ^  O  =  Us(S>°)  t3K  +  ~]  +  U(S,0)[K»  + 


nn '  '  '  nnO 

2Ksina  2sinacosa 


(3.53) 


From  (3.50),  (3.49),  (3.46)  and  (3.45),  we  have 


sma 


U  (s,0)  =  U  n  =  -U  (s,0)  -  U  (s,0)  + 

nnv  ’  nnO  ss^  sv  r^ 


mpivSrl  .  K,t 


•  2  , 
sina  +  sin  a] 


(3.54) 
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In  the  above  equations,  U(s,0)  is  the  solution  of  the  integral 
equation  (3.11),  and  hence  Us(s,0)  and  Uss(s,0)  are  calculated  nuneri- 
cally.  Having  thus  obtained  V„,  D#,  and  Un„  fron  equations  (3.51) 
through  (3.54),  we  can  now  calculate  U(s,n)  and  V(s,n)  fro*  equations 

(3.35)  and  (3.36). 
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CHAPTER  IV 

THICK  BOUNDARY  LAYER  AND  NEAR  WAKE 

This  chapter  describes  the  method  of  calculating  the  development 
of  the  thick  turbulent  boundary  layer  over  the  body  and  wake.  The  com¬ 
plete  Reynolds  equations  for  the  mean  flow  were  first  examined,  and 
suitable  approximations  were  made  to  simplify  them  in  [4].  The 
resulting  equations  are  elliptic  and  contain  a  number  of  additional 
terms  in  comparison  with  the  usual  thin  boundary- layer  equations. 

1 .  Governing  Equations 

The  momentum  and  continuity  equations  for  a  thick,  axisymmetric 
boundary  layer,  given  in  [  5]»  are 


h^  3x  3y  ph^  3x  rh^  3y  1  p  J 


L|V  +  v3V.Lu2  +  13£=  o 

hx  3x  V  3y  hj  U  p  3y  U 


(4.1) 


(4.2) 


h (Ur)  *  h  (rhiv)  =  0 


(4.3) 


where  U  and  V  are  the  mean  velocity  components  in  the  boundary  layer 
along  tangential  and  normal  directions  of  the  body,  p  is  the  pressure 


at  a  point,  h^  =  1  +  Ky,  is  the  metric  coefficient  in  the  longitudinal 
_  £1] 

direction,  x=  -p  uv  +  u  is  the  shear  stress  and  p  is  the 
dynamic  viscosity  of  the  fluid.  See  fig  [*.l]. 
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Equation  (4.2),  in  which  viscous  and  turbulence  terms  are  neglected, 
implies  that  the  viscous  and  turbulent  stresses  acting  in  the  normal  di¬ 
rection  are  negligible  in  comparison  with  the  pressures  and  the  variation 
of  static  pressure  across  the  boundary  layer  is  associated  primarily  with 
the  curvature  of  the  mean-flow  streamlines. 

2.  Turbulence  Model 

The  adopted  turbulence  model  discussed  in  [3]  is  written  here  for 
convenience.  The  equation  for  the  conservation  of  turbulent  kinetic 
energy  was  transformed  into  one  for  the  Reynolds  stress  t : 


1 /■  U  3t  y  3 tv  _  r  3U  t/ji\  . 

is;  (k:  »  * v  37)  -  T(s7  -  “»  * 


*'»[{«  i 


,  jn 
+  i-77i  =  ° 


(4.4) 


where  a^^  (=  - j)  is  a  constant  (=0.15),  G(y/6)  is  a  diffusion  function 

q 

and  t(y/6)  is  a  length-scale  function  identified  with  the  usual  mixing 


length . 


The  equation  for  ' £'  that  was  used  is  given  by 


aeeff 


(4.5) 


3x  =  <e  ~  eeff^106> 
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where  i  is  the  length  scale  with  the  usual  rate  of  strain  — . 
l  is  the  length  scale  with  the  extra  rate  of  strain  e  and  a  is  a 
constant  of  the  order  10,  eg^  is  the  effective  rate  of  strain,  106  re¬ 
presents  the  lag  length  over  which  the  boundary  layer  responds 
to  a  change  in  e. 

3.  Method  of  Solution 

The  axisymmetric  laminar  boundary  layer  from  the  nose  to  the  natur¬ 
al  point  of  transition  was  calculated  by  usingThwaite's  method.  From 
the  natural  point  of  transition  onward,  equations  (4.1),  (4.2)  and  (4.3) 
were  used  to  calculate  the  thick  turbulent  boundary  layer. 

Equations  (4.1),  (4.2)  and  (4.3)  contain  four  unknowns,  namely 
U,V,p  and  uv.  Even  when  a  turbulent  model  is  introduced  for  uv,  the 
resulting  set  of  equations  cannot  be  solved  by  a  marching  technique. 

This  is  due  to  the  presence  of  the  pressure  as  an  unknown. 

Equations  (4.1)  and  (4.3)  are  first  solved  by  a  marching  technique 
for  U  and  V  for  an  assumed  pressure  field  p(x,y).  The  y  momentum  equa¬ 
tion  (4.2)  is  then  usea  to  update  the  pressure  field.  In  this  itera¬ 
tion  scheme, the  solution  for  the  pressure  p  and  of  the  y  momentum  equa¬ 
tion  lags  the  velocity  field  by  one  iteration.  The  solution  is  ob¬ 
tained  through  step-by-step  integration  by  marching  downstream  from 
some  initial  station  where  the  velocity  and  shear- stress  profiles  are 

3U 

prescribed.  In  the  wake  calculation,  the  conditions  x  =  0  and  —  =  0 
at  the  wake  centerline  are  used. 
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CHAPTER  V 

INTERACTION- ITERATION  PROCEDURE 


To  begin  with,  the  velocity  distribution  on  the  body  is  calculated 
from  the  integral  equation  of  the  first  kind  for  the  unknown  velocity 
for  irrotational  flow  outside  the  body.  Using  the  above- calculated 
pressure  distribution,  the  laminar  boundary  layer  is  calculated  up  to 
the  natural  point  of  transition,  and,  from  then  onwards,  the  turbulent 
boundary  layer  is  calculated.  In  continuing  the  calculations  into  the 

t 

au 

wake,  the  conditions  —  =  0  and  t  =  0  are  applied  on  the  wake  centerline. 
The  first  time  the  boundary  layer  and  wake  are  calculated,  the  matching 
boundary  between  the  boundary  layer  and  the  external  irrotational  flow 
is  fixed  once  and  for  all  at  1.256  from  the  wall.  This  ensures  in  this 
particular  case  that  the  fixed  matching  boundary  always  lies  in  the  ex¬ 
ternal  irrotational  flow.  The  stream  function  at  the  matching  boundary 
is  calculated  and  is  used  to  calculate  the  normal  velocity  at  the  matching 
boundary.  This,  in  turn,  is  used  as  a  Neumann  boundary  condition  for 
the  external  potential  flow.  This  solution  is  then  employed  to  update 
the  pressures  along  the  boundary,  and  the  turbulent  boundary  layer  is 
recalculated.  This  iteration  procedure  is  carried  out  until  the  re¬ 
quired  convergence  is  achieved  in  the  wall  pressure  vtribution.  Each 
time, the  turbulent  boundary  layer  is  calculated  from  the  point  of  natur¬ 
al  transition,  where  velocity  and  shear  stress  profiles  are  given. 
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At  this  point  it  is  important  to  note  that,  due  to  the  coordinate 
transformation  at  the  first  station  of  the  wake  calculation,  the  results 
obtained  in  the  wake  are  incorrect  and  so  for  the  purposes  of  calculating 
the  drag,  a  different  control  volume  is  used,  as  is  described  in  Chapter 
VI,  instead  of  the  far-wake  momentum  thickness. 

In  the  case  of  a  sink  behind  the  tail,  the  external  irrotational 
flow  equations  were  modified  in  Chapter  III. 


'■‘•*”WSSSP^ 
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CHAPTER  VI 
CALCULATION  OF  DRAG 

The  purpose  is  to  compare  the  computed  values  of  the  thrust  de¬ 
duction  from  potential  flow  and  viscous  flow  methods.  Since  the 
drag  experienced  by  a  body  in  a  uniform  irrotational  flow  is  zero, 
the  calculated  drag  of  the  body  with  an  external  sink  equals  the  dif¬ 
ference  in  drags  ie.  the  "thrust  deduction".  The  calculation  of 
drag  in  potential  flow  will  be  discussed  in  Section  1.  The  difference 
in  drags  computed  with  and  without  a  sink  in  viscous  flow  gives  the 
value  of  the  thrust  deduction  in  viscous  flow.  This  procedure  will  be 
described  in  Section  2. 

1.  Calculation  of  Drag  in  Potential  Flow 

As  mentioned  earlier  [Ch.  l],a  calculation  procedure  for  drag 
in  potential  flow  that  is  accurate  by  at  least  an  order  of  magnitude 
higher  than  that  given  by  the  H.-S.  method  [  2]  is  required.  Such  a 
method  was  derived  and  programmed  by  Landweber,  L. ,  and  Choi,  D.H., 
and  will  be  presented  here  briefly. 

An  integral  equation  of  the  second  kind  in  terms  of  the  unknown 
vortex  distribution  (tangertial  velocity)  is  solved. 

Consider  the  case  of  the  body  stationary  in  a  uniform  stream  of 
velocity  u  =  1  in  the  negative  x-direction.  Assume  that  the  fluid 
within  the  body  is  at  rest,  and  that  the  disturbance  velocity  field 
is  generated  by  a  vortex  sheet  of  strength  y  on  the  body  surface. 


sy«*M 
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Let  u(s)  denote  the  tangential  velocity  at  the  body  surface  and 
y(s)  the  magnitude  of  the  vector  f,  where  s  is  arc  length  along  a 
meridian  section,  measured  from  the  tail.  See  fig  [6.1].  Then 

u(s)  =  -y(s)  (6.1) 


The  condition  that  the  tangential  velocity  on  the  interior  side 
of  the  body  surface  S  be  zero  yields  the  integral  equation 


1/2  y(s)  + 


1_ 

4tr 


^■yvv 


0  0 


r  dt  d6 


=  cos  a 


(6.2) 


PQ 


in  which  the  first  term  is  contributed  by  the  local  vortex  element, 
the  second  is  that  of  the  entire  vortex  sheet,  given  by  the Biot-Savart 
law,  and  the  third  is  due  to  the  uniform  stream.  Here  P  denotes  a 
fixed  point  at  (s,6p  =  0),  and  Q(t,8)  a  variable  point  of  integration 
in  the  arc  length  t  and  polar  angle  8.  The  total  length  of  arc  along 
a  meridian  is  i.  The  notations  r^  =  r(t) ,  rg  =  r(s) ,  rpQ  for  the 
position  vector  of  P  relative  to  Q  are  used.  The  magnitude  of  rp^  is 
given  by 


rPO  =  V  +  Tt  +  rs  -  2rsrt  C0S6] 


1/2 


tRst2  "  4  rsrtcos2e/2J1/2 


where  R 


st 


■  (xt  -  V  ♦  <rt  *  rs5 


and  a  is  the  angle  made  by  the  tangent  with  the  x  axis. 


r- 
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The  most  evident  problem  faced  in  solving  the  integral  equation 
is  the  numerical  treatment  of  the  singularity  of  the  kernel.  The  pro¬ 
duct  of  the  strength  of  the  vortex  sheet  and  the  secant  of  the  slope 
angle  was  chosen  as  an  unknown  variable  instead  of  vortex  distribution 
itself,  since  this  product  varies  rather  mildly  over  the  entire  body, 
whereas  the  secant  of  the  slope  angle,  which  appears  in  the  resulting 
integral  equation,  goes  to  infinity  at  both  ends  of  the  body.  Yet 
the  integrand  has  a  logarithmic  singularity  when  the  moving  point 
and  the  fixed  point  coincide.  This  can  be  removed  by  adding  an  iden¬ 
tity  to  the  integral  equation. 

After  evaluating  the  scalar  triple  product  and  putting 
r(s>  =  y(s)  sec  os,  eqn  (6.2)  becomes 


r(S)  =  2  "  4*  o  o  “3"  *  {2rsrt  "  £f,(Vxs)+2fs]  cosel  dedt 

rsrPQ 


(6.3) 


The  above  equation,  after  integrating  with  respect  to  e  and  eli¬ 
minating  the  singularity,  has  an  integrand  that  remains  finite  at  the 
body  ends,  and  is  of  the  form  of  a  product  of  a  slowly  varying  un¬ 
known  function  and  a  rapidly  varying  kernel,  involving  elliptic  inte¬ 
grals. 

To  solve  the  above  integral  equation,  the  entire  body  is  divided 
by  a  certain  number  of  evenly  spaced  sections  along  the  arc  length 
at  which  the  solution  was  sought.  In  order  to  achieve  higher  accuracy, 
however,  it  was  necessary  to  subdivide  each  section  further.  At  these 
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additional  node  points,  the  solution  is  not  sought  directly  but  ob¬ 
tained  from  the  four-point  Lagrange  interpolation  formula  by  taking 
advantage  of  the  fact  that  the  unknown  is  a  slowly  varying  function. 
Now  the  integral  equation  is  discretized  with  the  use  of  a  Simpson's 
rule  and  the  resulting  simultaneous  equations  are  solved  by  the  Gauss- 
Siedel  interation  procedure. 

After  the  calculation  of  T(s),  is  calculated  from  the  formula 


J 

d  0 


r2( 


4f+f 


3/2 

r)  f’ds 


(6.4) 


2 

where  d  is  the  maximum  diameter  and  r  =  f(x)  is  the  equation  of  the 
body. 

At  each  stage,  the  method  was  applied  to  a  spheroid  of  length- 
diameter  ratio  Jl/d  =  6.0,  and  the  resulting  velocity  distribution  was 
compared  with  values  from  the  known  analytical  solution.  The  pres¬ 
sure  drag  coefficient  was  computed  for  only  half  the  spheroid  for 
which  an  exact  analytical  formula  is  known.  The  computed  value  was 
CQ  =  -  0.00857572;  the  exact  value  is  -0.00857578  and  this  accuracy 
is  considerably  better  than  required  for  the  thrust-deduction  problem. 

The  above  program  was  modified  for  an  axial  downstream  sink  in 
a  uniform  flow  field  and  the  drag  on  the  body  was  computed  for  the 
EPH  form  with  a  non-dimensional  sink  strength  of  M'  =  -0.004  (Ch.  2) . 
The  computed  drag  coefficient  is  0.00392. 


36 


2.  Calculation  of  Viscous  Drag 

Due  to  the  coordinate  transformation  at  the  first  station  of  the 
wake  calculation,  the  results  obtained  in  the  wake  are  not  very  accur¬ 
ate  and  so  for  the  purposes  of  calculating  the  drag,  a  special  control 
volume  is  used.  This  control  volume  is  shown  in  fig  [6.2].  The  body 
is  considered  to  be  at  the  center  of  a  circular  channel  of  large  ra¬ 
dius,  with  the  flow  considered  to  be  inviscid  in  the  neighborhood  of 
the  channel  wall.  It  is  convenient  to  select  as  a  boundary  of  the 
control  surface  a  transverse  section  perpendicular  to  the  x-axis,  AB, 
far  ahead  of  the  body  and  a  transverse  ring  made  by  CD,  meeting  another 
ring  made  by  DE  and  enclosing  the  remaining  tailpiece  EF  of  the  body. 

E  is  the  last  point  of  calculation  on  the  body  surface  and  DE 
is  the  normal  to  the  body  meeting  EBLW  at  the  point  D.  In  the  boundary 
layer,  the  boundary  of  the  control  volume  is  chosen  to  be  along  the  normal 
to  the  body  because  the  boundary  layer  calculations  are  done  in  a 
tangential  and  normal  coordinate  system,  and  all  the  quantities  are 
known  along  normals  to  the  body. 

The  momentum  theorem  gives  the  expression  for  the  drag  D  of  a  body 
in  a  steady  flow, 

D  =  -  /  [pfc  +  pu(y*n)]  dS  (6.5) 

S 

as  an  integral  over  a  control  surface  S  surrounding  the  body,  where 
p  is  the  pressure,  n  is  the  outwardly  directed  normal,  a,  m  and 
n  are  its  direction  cosines,  V  is  the  velocity  vector,  ug  the  s 
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component  of  the  velocity  vector,  u  the  x  component  of  the  velocity 
at  points  of  the  control  surface  and  p  the  density  of  the  fluid,  assumed 
to  be  constant.  Here  the  flow  incident  upon  the  body  and  the  direc¬ 
tion  of  the  drag  force  are  taken  in  the  positive  x-direction. 

In  applying  the  momentum  equation  (6.5)  to  the  selected  control 
volume,  we  see  that  there  is  no  contribution  to  the  integral  from  the 
boundary  surface  of  the  channel.  We  obtain,  then. 


A  procedure  for  obtaining  an  expression  for  the  drag  as  an  inte¬ 
gral  extending  only  across  the  wake  was  suggested  by  Betz  [6] .  Let 
us  suppose  that  the  potential  flow  outside  the  boundary  layer  and  wake 
is  continued  analytically  into  these  regions.  This  defines  a  potential- 
flow  field  with  velocity  components  u^,v^  *nd  pressure  p^,  which  dif¬ 
fers  from  the  original  flow  field  only  within  the  boundary  layer  and 
wake.  The  Bernoulli  equation  may  be  applied  to  obtain 

pi  =  po  +  I  p(u»  "  uis2)  (6>7) 

The  increase  in  flux  of  fluid  across  the  plane  FG  is  attributable 
to  sources  bounded  by  the  control  surface  of  total  strength 

tuls  '  Us)dS2 

b2 


m 


(6.8) 
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which,  gives  rise  to  a  negative  drag,  (subscript  s  denotes  tangential 
quantities) 


Ds  =  -Pj  Ujuls  .  us)dS2  (6.9) 

s  2 

Also,  we  have  for  the  flow  system  generating  the  potential  flow  u  , 

Vj  from  the  momentum  equation  (6.5), 


where  rfa  is  the  body  ordinate  at  the  last  point  of  calculation  and 
Pb  is  the  pressure  on  the  body  at  this  point,  assumed  to  be  constant 
near  this  region.  Subtracting  equation  (6.10)  from  (6.6)  gives 

D  =  Ds  +  /  [(Pj-p)*  +  pOY»1s  "  uuj]  dS2  +  **b2(plb-pb)  (6.11) 

s2 

Since  outside  S  ,  p1  =  p  u.  =  u  and  v..  =  v. 

JL  y  X  X 

Substituting  for  p^  and  Dg  from  (6.7)  and  (6.9)  into  (6.11)  gives 

D=  4/  {(PQ-P)  +  ^p[U2-uls2]}dS2 

s2 

+  p|  fusCU»u>  '  uis(U«ul)]dS2  +  "b^lb'V  (6‘12) 

S2 


an  expression  for  the  drag  which  involves  only  quantities  in  the 
boundary  layer. 

Now  a  sink  of  strength  M  is  introduced  in  the  wake  on  the  axis 
VpX  and  v^r  are  the  velocity  components  along  x  and  r  of  the  net  velo¬ 
city  due  to  the  sink  and  its  image  system  as  calculated  from  potential 
theory  [Ch.  3,  Sec  3].  Let  the  quantities  computed  with  the  sink  in 
the  wake  be 


u  *  u\  V  ♦  v*.  Vu|,  p  ♦  p*.  plb  *  p*b,  pb  *  p* 

U1  *  vpx'  V1  *  V’  “ls  *  VPS 


where  *  indicates  the  small  increments. 

The  drag  in  the  presence  of  a  sink  can  now  be  written  as 


up0-p-p*)  *  |  ti£  -  (uls .  vps)2]ids2 


+  pf  [(u  +  u*)(U-u-u*)-(u.  +v  )(U-u..-v  )]dS 
1  s  s  1  v  ls  ps  ®  1  px  J  2 


+  ,rb  (pib  *  ^b  -  Pt  -Pb’ 


(6.13) 


Subtracting  (6.12)  from  (6.13)  and  neglecting  higher  order  terms,  we 
obtain  for  the  difference  in  drags 


D'  -  D  "  -*/  (P*  +  PuisVps)dS2  +  Urb  (P!b  '  Pb} 


L  [U(us  • 


+  (u.v  -  uu*)]  dS„ 
1  ps  s'J  2 


(6.14) 
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CHAPTER  VII 
DISCUSSION  OF  RESULTS 


The  drag  coefficient  C^  of  the  body  based  on  wetted  surface  area, 
as  calculated  from  equation  (6.12),  is  0.00212.  This  estimate  is  com¬ 
pared  with  the  modified  Squire-Young  formula  of  Ref.  [7],  equation  (10), 


A„  U  1/8  [7(A7  +  2)  +  « 
2  _2  (-«)  2 
S  lU  J 


where  C  =  - D  is  the  drag  coefficient  based  on  wetted  surface  area 

U  l/2pujs 

S,  A^  is  the  displacement  thickness,  A2  is  the  momentum  thickness,  U^ 
the  reference  velocity,  U  the  velocity  at  the  edge  of  the  boundary  lay- 

v 

er.  All  the  above  quantities  are  to  be  evaluated  at  the  tail  of  the  body 
Since,  in  our  case,  the  boundary- layer  calculation  stops  before  the 
tail  of  the  body,  the  values  calculated  at  the  last  station  on  the  body 
are  used.  This  should  give  only  a  small  error  since  the  last  point  is 
very  close  to  the  tail.  Applying  the  computed  values 


Aj  =  0.133832  ft2 

A  =  0.090478  ft2 

A1  Ue 

-r±  =  1.47917  and  ~  =  0.93226 

A2  U« 

S  =  59.6588  ft2 

we  obtain 


CD  =  0.00239 


This  value  exceeds  that  from  (6.12)  by  12.8  percent.  This  dis¬ 
crepancy  may  be  due  to  the  small  change  in  momentum  thickness  from  the 
last  station  to  the  tail,  which,  according  to  Lee  [3],  Fig.  26. g,  would 
decrease  the  above  value  by  about  5  percent,  leaving  a  discrepancy  of 
about  8.5  percent. 

The  calculated  value  of  the  thrust  deduction  from  the  potential- 
flow  method  is  0.00021,  where  as  the  value  obtained  from  the  boundary- 
layer  potential-flow  interaction  method  is  0.00043.  This  calculation 
can  be  improved  by  calculating  the  effect  due  to  the  sink  more  accurately; 
i.e.,  by  using  a  more  accurate  set  of  inputs  for  Green's  third  formula. 

The  boundary- layer  and  wake  equations  need  to  be  modified  for  a  case 
where  a  singularity  is  present  in  the  flow  field.  Also  the  procedure 
for  calculating  the  boundary  layer  and  wake  needs  to  be  improved  so  as 
to  obtain  a  smooth  transition  from  the  boundary  layer  to  the  wake. 
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APPENDIX 

To  prove  that  the  term  d<j>o  +  is  an  exact  differential. 

We  have 


HQ 


3*o 


%  d*o  +  *d  d*o  =  %  C1T  dx  +  ir  dr)  +  *d(  aT  dx  +  "aT  dr) 


3*o  H>  84,o  **o 

l(*dF*  *i  » r> dx  *  <%  a-  *  <a  ir1  dI 


a) 


The  condition  that  this  be  an  exact  differential  is  that 


3r  ^d  3X  3x 


9  3$ 


8*, 


o 

fil;  -  +  41  j  V  ) 

9*  vvd  3r  d  9r 


(2) 


or 


^dr  ^ox  +  ^d  ^oxr  +  ^dr  ^ox  +  ^d^oxr  ^dx  ^or  +  ^d  ^oxr  + 


*dx  tv  +  *d  ixr 


or 


^dr  ^ox  +  ^ir  ^ox  ~  ^dx^or  +  ^dx^or 


(3) 


We  have  the  following  relations  between  the  potential  and  stream 
functions  in  cylindrical  coordinates: 

3(frd  _  .  3^°  - 

r  1?  3x  ’  r  Hr  "  ax  (4) 

_  3lHi  3+o  3*o 

r  lx  ar  ’  r  lx  9r 


'.'-.imMww 


5^ HEgg^ipi 
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Hence,  we  obtain  d>,  *  =  (r  6  ,  )  f— »D  )  =  cb .  <p  and  similarly 

’  Ydr  Yox  v  rdx  r  ror'  Ydx  or  1 

^dr  ^ox  =  ^dx  ^or*  substituting  these  relations  in  (4)  we  realize 
that  condition  (2)  is  satisfied  and  hence  ^  d^  +  ^  d^  is  an  exact 


differential. 
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